Supersolidity in the triangular lattice spin-1/2 XXZ model: A variational perspective 
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We study the spin-1/2 XXZ model on the triangular lattice with a nearest neighbor antiferromag- 
netic Ising coupling J z >0 and unfrustrated (J± < 0) or frustrated (J± > 0) kinetic terms in zero 
magnetic field. Incorporating long-range Jastrow correlations over a mean field spin state, we obtain 
the variational phase diagram of this model on large lattices for arbitrary J z and either sign of J±. 
For J± < 0, we find a vo supersolid for J z /\ J±\ > 4.7, in excellent agreement with quantum 
Monte Carlo data. For J± >0, a distinct %/3x\/3 supersolid is found to emerge for J z /J± > 1. Both 
supersolids exhibit a spontaneous density deviation from half-filling. At J z /J±=oo, the crystalline 
order parameters of these two supersolids are nearly identical, consistent with exact results. 



PACS numbers: 75.10.Jm 05.30.Jp 71.27.+a 

Introduction. — Understanding how quantum effects 
select a unique ground state from an exponentially large 
number of classically degenerate configurations is a prob- 
lem that is common to geometrically frustrated magnets, 
Mott insulators, and quantum Hall systems. The sim- 
plest example of a frustrated classical magnet is the two- 
dimensional (2D) classical Ising antiferromagnet on the 
triangular lattice, which has an extensive ground state 
entropy, with an entropy per spin ps 0.323, and critical 
spin correlations [lj. The application of an in- 
finitesimal transverse magnetic field is known to lift this 
ground state degeneracy and favor a long-range v3 x vo 
ordered state [2j where the sublattice magnetizations on 
the three sublattices take on the form (0, m z , —m z ). An- 
other quantum variant is the triangular lattice spin-1/2 
XXZ model, described by the Hamiltonian 



H 



Jx(S*Sl 



(1) 



where J z > and (ij) refer to nearest neighbour links of 
the triangular lattice: interest in this model stems from 
an early suggestion |3J] that this model might support, for 
J_L > 0, a quantum spin liquid ground state for J z /J± >• 
1. (Henceforth, we will set |Jj_| = 1.) This model may 
also be viewed as a hard-core boson model with a nearest 
neighbor boson repulsion and a boson hopping amplitude 
which is unfrustrated (Jj_ < 0) or frustrated (Jj_ > 0). 
Interest in such models has been revived following reports 
[H of supersolidity in 4 He which remain to be understood. 

In this paper, using a variational wave function which 
incorporates long range Jastrow correlations, together 
with recently developed Monte Carlo optimization meth- 
ods we obtain the variational phase diagram of this 
model on large lattice sizes for either sign of J±. The 
main significance of our work, apart from the fact that it 
represents substantial progress in the study of the XXZ 
model, is that it illustrates the power of the variational 
wave function approach in studying strongly correlated 
phases and quantum phase transitions. 



We begin by summarizing our key results in the context 
of recent important work on this model. 

For J±<0: (i) We incorporate Jastrow correlations 
above a uniform superfiuid wave function and find, using 
the Binder cumulant, a phase transition from a corre- 
lated superfiuid state into a supersolid with long range 
\/3 x v3 crystal order for J z >4.7; this supersolid order 
persists even for J z — ¥ oo. This is in good agreement 
with quantum Monte Carlo (QMC) simulation results 
which found a superfluid-supersolid transition at 
J z w 4.5. (ii) We find that the Jastrow factors exhibit a 
long range 1/r tail in the superfiuid as well as the super- 
solid which is consistent with a linearly dispersing sound 
mode, (iii) The crystal order in the supersolid emerges 
as a spontaneously broken symmetry in our wave func- 
tion as a consequence of Jastrow correlations, but the 
Jastrow factors do not completely fix the crystal pattern. 
Incorporating a variational parameter which directly cou- 
ples to the difference between the (2m z , — m z , —m z ) and 
(0,m 2 ,— m z ) candidate crystal orders, we find that the 
former has a slightly lower energy, consistent with QMC 
results |H [TO, El as well as a recent variational study 
[l2T | of the quantum dimer model at J z = oo. (iv) We find 
that a supersolid wave function which permits for a small 
density deviation away from half-filling has lower energy. 

For J± > 0: (i) We incorporate Jastrow correlations on 
top of a state with coplanar 3-sublattice magnetic order. 
Using this wave function, we find that there is 120°-xy or- 
der for J z < 1 whereas there is supersolid order for J z > 1, 
with (2m z , — m z , —m z ) crystal order, in agreement with a 
recent suggestion based on analyzing the J z = oo limit Q 
and density matrix renormalization group (DMRG) re- 
sults on small system sizes [3]. (ii) The formation of the 
supersolid state is accompanied by a spontaneous density 
deviation away from half-filling. Equivalently, the XXZ 
model in zero field develops a spontaneous S z magneti- 
zation; this magnetization grows with increasing J z . (iii) 
The 120°-xy state as well as the supersolid exhibit 1/r 
decay of Jastrow factors consistent with a linearly dis- 
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persing excitation mode, (iv) For J z — oo, recent work 
has shown,through the construction of a unitary trans- 
formation [6( as well as by a high temperature expansion 
Q, that the diagonal spin correlations are identical for 
either sign of Jj_ hi this limit, which enables one to con- 
clude that the frustrated XXZ model must also exhibit 
\/3x -y/3 crystal order, and is likely to be a supersolid at 
J z = oo 0,0j- At J z = oo, the crystal orders in our 'frus- 
trated supersolid' and the 'unfrustrated supersolid' are 
nearly identical and close to QMC results for J± < 0, in 
agreement with this result, (v) In contrast to the result 
of Ref . @ , we find that the off-diagonal order parameter 
for J± > is much smaller than that for J± < 0. 

Variational Wave function. — We begin with a varia- 
tional wave function of the form 

|lM - Je exp(]>>(r)S.fSf +r )7>G|^MF), (2) 

where |"0mf) is a variational mean-field wave function 
for spin-f and spin- J, fermions (see below), Vg is the 
Gutzwiller projection operator which restricts the Hilbert 
space to one fermion per site thus taking us from a 
fermionic state to a spin wave function, J(r) are vari- 
ational Jastrow parameters which build correlations be- 
tween spins separated by r while respecting lattice sym- 
metries, and Jq is a global variational parameter, elabo- 
rated upon later, which fully fixes the crystalline order. 

We choose I^mf) to be the ground state of a mean 
field Hamiltonian H MF = - ^4^-^ ft™ r -4 

where Si = \c\ a a a pCjp. The first term in 7?mf represents 
nearest neighbor hopping of fermions, with amplitudes 
i™ r , while the second term represents a site-dependent 
magnetic field whose direction and magnitude can be 
chosen to yield various magnetically ordered states. Note 
that the wave function has the freedom to decribe a mag- 
netically disordered spin liquid state if ft.™ 1 = 0. 

For the unfrustrated model (Jj_ = —1) we set t™ = 0, 

and select two different configurations of h™ . One choice 
is to set ftj ar = x, which leads to a mean field ferromagnet 
with the moment polarized along the S x direction. This 
mean field state has no variational parameters; we refer to 
the spin wave function, after Gutzwiller projecting this 
state and inclusion of correlation factors, as |>I'y M_x ). 
The second choice, which we motivate below, promotes 
3-sublattice magnetic order, with ft,™ 1 ' = (1j — 1/2, — 1/2) 
on the 3 sublattices, ft™. 1 = h x , and ft™ r = 0, where h x is 
a variational parameter. We refer to the correlated spin 
state obtained from this mean field state as |*™" xz ). 

For the frustrated model, we choose the hopping am- 
plitudes such that |i™ r | = 1 but with signs picked such 
that upward (downward) pointing triangular plaquettes 
enclose flux 7r (zero). Upon Gutzwiller projection, this 
state describes a spin liquid which preserves all lattice 
symmetries. It is equivalent to a pairing wave function 
which is known, from earlier studies Q to be a good start- 




FIG. 1: (Color online) (a) Order parameters, extrapolated 
to the thermodynamic limit, for J± = — 1. (b) Jastrow pa- 
rameters showing 1/r behavior at long distance (curves for 
J z — 4, oo have been offset downward for clarity). 

ing point to understand the triangular Heisenberg anti- 
ferromagnet. Motivated by the classical ground state, 
we pick ft™ 1 to describe 3-sublattice coplanar magnetic 
order. The spin wave function obtained from choos- 
ing fcg* = (h,-h/2,-h/2), ft™ 1 ' = (0,W3/2,-W5/2), 
and ft,™ 1 = 0, with ft being a mean field variational pa- 
rameter, will be called xy ). The spin wave func- 
tion obtained from choosing ft™ r = (0,h x ,—h x ), ft™ r = 
(h z , —h z /2, —h z /2), and ft™ 1 = 0, with mean field varia- 
tional parameters h x ,h z , will be called \^y F ~ xz ). 

In order to obtain the best trial wave function, we opti- 
mize the variational parameters using Monte Carlo sam- 
pling of the spin configurations combined with recently 
developed techniques [5| . We begin by discussing our re- 
sults for the unfrustrated model, for which we present 
comparisons with existing QMC results. We follow it 
up with results on the frustrated model which is less well 
studied since the sign problem does not permit controlled 
QMC simulations. Finally, we consider J z = oo, and 
demonstrate consistency with exact results in this limit. 

Ji < 0, Unfrustrated Case. — For J z -C 1, the ground 
state of the unfrustrated model is a quantum ferromagnet 
or, equivalently, a superfluid state of hard core bosons. 
With increasing J z , the ground state of the superfluid 
becomes increasingly correlated, and QMC simulations 
[§-11 1 show that a supersolid ground state, with coexist- 
ing superfluid and charge order, develops for J 2 >4.5. If 
we set Jo — 1 and optimize | , I'y M ~ x ), we find that the 
ferromagnetic order parameter m 2 ^ = |(S ,:E )| 2 is nonzero 
for all values of J z , as seen from Fig. [lja), consistent 
with superfluidity persisting for arbitrarily large inter- 
actions. For all interaction strengths, we also find, as 
seen from Fig. [IJb), that the optimal Jastrow parame- 
ters J{v) decay as 1/r with distance. Such a power law 
decay implies a density structure factor scaling as ~ |q| 
at small momenta which is consistent with a linearly dis- 
persing Goldstone mode above a superfluid ground state 
hj. it is a nontrivial test of the unbiased variational 
optimization that we recover this behavior by optimizing 
over a very large set of variational Jastrow parameters. 
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FIG. 2: (Color online) (a) Binder cumulant crossing for m z 
displaying a superfluid-supersolid transition at ji crlt ' 1 « 4.7. 
(b) Probability distribution of the phase of order parameter 
for he — (dashed line) and at the optimal he>0 (solid line) . 



In order to probe for supersolidity, we compute the 
density structure factor (equivalently S z correlations). 
Although I*!/™ - *) is a translationally invariant state, 
we find that the density structure factor diverges at mo- 
menta ±(47r/3, 0) for large J z which shows that there is 
\/3 x \/3 crystalline ordering, that coexists with the su- 
perfluidity, for strong interactions. It is remarkable that 
the crystalline order emerges as a spontaneously broken 
symmetry induced by Jastrow correlations in our wave 
function. Using the crossing of the Binder cumulant on 
different system sizes as shown in Fig. [UJa), as well as 
from a calculation of the order parameters extrapolated 
to infinite system sizes shown in Fig.QJa), we locate the 
superfluid to supersolid transition point at J z ~ 4.7, 
which is in good agreement with QMC simulation results. 

While the crystal order is spontaneously generated in 
this wave function, a complete characterization of the 
crystalline order requires knowledge of the phase of the 
order parameter $ — S\ + S%e l27T ^ 3 



.£<z e i47r/3 = |$| e * 



where S\ B c refer to the total density on the A, B, C 
sublattices. Specifically the two candidate v3 x \/3 
crystal orders for this model have 9 = 2mt/ 6 with 
n = 0, • • • , 5, which corresponds to having a density 
order (2m z , —m z , —m z ) on the three sublattices, and 
9 = (2n + l)7r/6 with n = 0, • • • ,5, which implies 
(0, m z , —m z ) order. Since the Jastrow factors are insen- 
sitive to the phase of $, they cannot select between these 
candidates. In order to find the optimal supersolid, we 
therefore include a factor Jg = exp(hg cos 66*) in the wave 
function, with hg being a variational parameter. Opti- 
mizing this variational parameter, we find that hg > 
has the lowest energy indicating that the best supersolid 
state has (2m z , —m z , —m z ) crystal order. As shown in 
Fig. EJb) for J z — oo, the probability distribution P(8), 
of the phase 9, shows sharp peaks at the expected angles 
at the optimal hg > [14j , while the state with hg = is 
relatively featureless (the small peaks at 9 = (2n + l)7t/6 
appear to be a finite size artifact). Our results are con- 
sistent with earlier QMC simulation results [H El 
and recent variational calculations fl2| in the quantum 
dimer model limit at J z = oo. Motivated by these re- 



TABLE I: Exact diagonalization energy versus variational 
energy of |\f r y F-xy ) for a 6x6 system with J z < 1. 



J Z 


0.0 


0.5 


0.9 


1.0 


VMC 


-0.40950(7) 


-0.46992(5) 


-0.53075(5) 


-0.54816(9) 


ED 


-0.410957 


-0.475366 


-0.541988 


-0.560373 




FIG. 3: (Color online) (a) Energy difference AE = 



£* F_xy showing that \rj) 



AF-xz 
V 



) has lower energy for J z > 1 

AF-xy 



(b) Jastrow parameters versus distance for J z — 1 {\ip v 
and Jz = 3 (|V , v F_XZ ))- Solid line depicts 1/r behavior, 
for J z — 1 has been shifted upward for clarity. 
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suits, we study the properties of |4 , ™ _xz ) which imposes 
the above supersolid order, and find that the optimized 
wave function of this type exhibits a small density devia- 
tion away from half-filling. Equivalently, the XXZ model 
has a spontaneous S z magnetization in the supersolid. 

Ji > 0, Frustrated Case. — For the frustrated model, 
the sign problem prevents controlled QMC simulations on 
large lattices. There have been some classical and spin 
wave analyses for large spin 



15j . and exact diagonaliza- 



tion (ED) studies on small systems which do not ad- 
dress the issue of the thermodynamic limit for spin-1/2. 
The classical ground state with J z < 1 has 120° mag- 
netic ordering in the xy-plane. We therefore begin with 
the wavefunction |\[' V r F-xy ) as a variational candidate. 
With increasing J z , the system becomes more strongly 
correlated, developing strong short range S z correlations. 
Table I shows that the variational energy compares well 
with ED results [l6| on a 6 x 6 lattice for J z < 1. 

From the density structure factor and probability dis- 
tributions of the phase of <E>, defined earlier, we find that 
this wave function also displays a spontaneously ordered 
V3 x y/3 supersolid, with staggered charge and super- 
fluid orders, for J z > 1.2. The charge order is of the type 



(2m 2 



s ) on the three sublattices, identical to 



what we found in the unfrustrated model, and m z grows 
continuously beyond J z ~1.2, whereas the superfluid or- 
der is of the type (0, mj_, -mi) on the 3 sublattices. Mo- 
tivated by the proximity of this transition to the Heisen- 



berg model, we compare the two states, |^ 



AF- 
V 



AF-xy\ 
V 



and 



to see if the plane of the coplanar order rotates 
at this SU(2) symmetric point. As shown in Fig. 
the energy of \^y F ~ xz ) is indeed lower for J z >1. We find 
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FIG. 4: (Color online) (a) Diagonal and off-diagonal order 
parameters extrapolated to the thermodynamic limit showing 
the order parameter discontinuity at J z — 1 and the coexis- 
tence of both orders for J z > 1. (b) Order parameters for 
a 6 x 6 system compared with DMRG results from Ref. 0|. 
(c) Density distribution in |^y F_xz ) showing nonzero mag- 
netization with the onset of supersolid order, (d) Diagonal 
and off-diagonal order parameters in the limit J z —<x show- 
ing good agreement of the extrapolated m\ for J± > with 
QMC data on the unfrustrated model. 



that the Jastrow factors in | ^y F xy ) 



as well as |*y F ~ xz ) 



exhibit a long range l/r tail consistent with a linearly 
dispersing mode (see Fig. EJb)). 

Using the two different wave functions on the two sides 
of the Heisenberg point, we conclude that m z and m± are 
discontinuous at J z — l, and |'Fy F_xz ) describes a super- 
solid state as seen from Fig. 4(a). For a 6x6 system, our 
results are in excellent agreement with DMRG calcula- 
tions [?J (see Fig. 4(b)). The crystal order grows mono- 
tonically for J z > 1, while m± gets strongly suppressed. 
We find that |'J , y F_xz ) exhibits a small spontaneous den- 
sity deviation from half-filling, as shown in Fig. 4(c), 
which grows with J z to a maximum value 0.05%. 

J z = oo, the Quantum Dimer Limit. — At J z = oo, 
the Hilbert space of the XXZ model gets constrained 
to the minimally frustrated classical configurations of 
the triangular lattice Ising antiferromagnet Remark- 
ably, although the 'frustrated supersolid' wave function 
|\I/y F_xz ) describes a very different state from the 'unfrus- 
trated supersolid' | , F FM ~ X ), we find, upon projecting to 
the minimally frustrated Ising subspace, that the crystal 
order parameters of both states are nearly equal, and are 
in good agreement with the extrapolated QMC result for 
the unfrustrated case, as shown in Fig.HJd). We also find 
that the energies of these two completely different states 
are fairly close, within 6% of each other. This is consis- 
tent with recently uncovered exact results @, 0] at J z = oo 



which show that the crystal order parameters and ground 
state energies should be equal for either sign of J±. As 
seen from Fig. 4(d), m± extrapolates to a small value for 
either sign of J±, with m±( J± = 1) <Cto_l( J± = — 1). 

Summary. — In summary, we have used a variational 
approach to obtain the phase diagram of the frustrated 
and unfrustrated spin-1/2 XXZ models on the triangu- 
lar lattice for large lattice sizes and arbitrary J z > 0; 
these regimes are not accessible via other techniques. The 
phase diagram of the unfrustrated model is in good agree- 
ment with QMC data. The frustrated model exhibits a 
wide regime of supersolidity with no evidence of a spin 
liquid. Our results for J± < could be tested using po- 
larized dipolar bosonic molecules in a triangular lattice 
by tuning the ratio J±/ J z . For the frustrated model, pro- 



posed schemes to induce fictitious gauge fields [17| need 
to be adapted to produce a frustrated molecule hopping. 
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